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Collective behavior of active elements inspired by mass of biological organisms is addressed. Espe- 
cially, two topics are focused on among amazing behaviors performed by colony of ants. First, task 
allocation phenomena are treated from the viewpoint of proportion regulation of population between 
different states. Using a dynamical model consisting of elements and external "stock materials", 
adaptability against various disturbances is numerically studied. In addition, a dynamical model 
for a colony ants interacting via two kind of pheromones is studied, in which simulated ants, as a 
mass, are shown to make an efficient foraging flexibly varying the foraging tactics according to feed- 
ing schedules. Finally, experiments are performed with robots moving in virtual pheromone fields 
simulated by CG and CCD camera feedback system. Trail formation processes are demonstrated 
by this multi-robot system. 



I. INTRODUCTION 



Among intriguing phenomena exhibited by biological 
mass, it is remarkable that various functions or perfor- 
mances of the whole mass exceed simple superposition 
of each individual's ability. The mass of individuals we 
consider here is that of slime molds, (social) insects, fish, 
birds, etc. For example, a kind of cellular slime mold 
shows regulation of ratio between two different types of 
cell during its differentiation process 0, S S B • About 
the motion of the group, birds flock, fish school and 
mosquitoes swarm J|, |^ |3, [l3 • Social insects exhibit 
lots of amazing performances like colony formation, age 
polyethism, group foraging, etcTl^. It is worth noticing 
that these collective performances do not require special 
individual(s) which controls the behavior of the whole. 
Hereditarily homogeneous individuals achieve these col- 
lective behaviors by interacting each other through direct 
sight or chemical materials like c-AMP, pheromones, etc. 
How does each of them "feel" the situation of the whole 
and how does it determine its behavior? Here, we focus 
on the task allocation phenomena and the pheromone 
trail formation processes among various collective behav- 
iors described above. 

This article is organized as follows. In the next sec- 
tion, from the viewpoint of proportion regulation, multi- 
task allocation phenomena are treated with a dynamical 
model consisting of simple elements and external vari- 



ables. In §3, foraging behavior of modeled ants was 
treated in which ants leave and follow only the local 
information of pheromone, however they can fulfill the 
highly efficient foraging flexibly adapting to various feed- 
ing conditions. In §4, real robots experiments are per- 
formed. Using a virtual pheromone system (V-DEAR) 
which consists of LC projector and CCD camera, both 
task allocation process and trail formation process are 
examined by multi-robot system. Summary is given in 
S5. 



II. TASK ALLOCATION 

Division of labor {polyethism) performed by social in- 
sects or cell differentiation observed in the life cycle of 
Dictyostelium discoideum exhibit typical proportion reg- 
ulation phenomena. Namely, individuals (or cells) with 
the same inheritance divide into two or more groups and 
they play different roles. The proportion of populations 
of each group is regulated against external disturbance 
without special individual. Regulation means a kind of 
adaptability, i.e., for large disturbance like loss of some 
part of colony or change environment, individuals know 
the change of the situation about the whole and change 
their behavior. As a result of such autonomous process, 
the proportion of population between the groups is main- 
tained. Several questions naturally arise for these phe- 
nomena, i.e., how do they know about the whole and 
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FIG. 1: Examples of rule between multi-states, (a) sequential 
(b) intermediate, and (c) democratic. Arrows denote the path 
which the rule allows to change. Note that the paths between 
the states in rule (a) and (c) are all symmetric while it is not 
in (b). 



how do they determine their behavior? In this section, 
dynamical models for internal states of individuals and 
external variables of chemical materials are used to re- 
produce these phenomena and consider their mechanism 
through numerical simulations. 

As a previous work, we suggest a dynamical model for 
proportion regulation phenomena for two states. The 
model consists of N identical elements with two internal 
variables, i.e., activator Ui{t) and inhibitor couple 
globally through their average quantities 0. The dy- 
namics of the internal variables Ui{t) and Vi{t) is 

Ui = aui — hvi — u\ ^ Ki{u — Ui)^ (1) 
Vi = cui - dvi ^ K2{v - Vi), (2) 

where ~ denotes averaging over population, i.e. , u = 
Uj/N and v = ^ ^j/^- and d are non- negative 

parameters. Ki and K2 are the non- negative suscepti- 
bilities of each variables in the context of fast diffusion 
limit. The model has following properties under suit- 
able parameter conditions a — d < 0, ad — he < and 
Ki = 0: (i) The population "differentiates" into two 
states for K2 > i^2c- (ii) The population ratio between 
the two states is regulated within some range which is ob- 
tained analytically, (iii) The system maintains the ratio 
against large disturbance by changing the state of indi- 
viduals, (iv) State of each element depends only on its 
initial condition. The model, however, is only available 
for two states allocation because of its simplicity. In order 
to treat multi-states case, we adopt another dynamical 
model in the next subsection. 



A. Model 

Unlike the two states case, multi-states case requires 
detail setting up the problem, i.e., not only the ratio 
between the states but also the rule between the states 




FIG. 2: 

states. 



Typical shape of potential function U with three 



are needed to be assigned. Here, the rule is considered 
to be determined by the connectivity between the states. 
Several examples of rule are shown in Fig^ We assume 
the sequential and symmetric rule (FigH^a)) here, and set 
up the situation as follows: identical individuals, M- 
states with proportion of population for j-th state with 
0<pj <1 and = 1. M = 3 and pj = {0.5, 0.3, 0.2} 
are used as typical setting up. 

Our model consists of A^ identical elements and M 
"stock materials". The elements represent the individ- 
uals and have internal variable Ui{t), i = 1,...,A'. The 
stock materials correspond the states and their quantities 
are expressed by external variable Wj{t) , i = 1,...,M. 
Each state is expressed by a valley of potential function 

U{u) = ^^{aj - Wj{t))Uj{u) + aoUo{u), (3) 



g-/3(«-S3) 

Uo{u) = w^, 



(4) 
(5) 



where sj and aj denote the positions and depths of j- 
th valley, respectively (Fig|2|). P gives an average depth 
of the valleys. These quantities are determined by the 
environment. Uq is used for boundedness of the internal 
variables of individuals. Wj{t)^ the quantity of j-th stock 
materials , changes effective depth of j-th potential and 
obeys following equation 



■Ujiui)), 



(6) 



where the first term of the R.H.S represents decay process 
and the second term is production process by the individ- 
uals staying in the corresponding state, tj is characteris- 
tic time of j-th material and a is a production rate. The 
dynamics of Ui{t)^ the internal variable of i-th element, 
is 



Ui = 



dU{u) 
du 



(7) 



where 77 is a white uniform noise with amphtude v. In- 
dividuals are considered to be subjected under the time 
dependent potential U {u) and the noise rj. 
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FIG. 3: State of j-th element uj (open rectangle) and the 
shape of the potential at t = (a) and t = 100 (b). The 
potential is modified by stock materials wj. To display the 



modification by Wj, initial shape of the potential (^ 



0) 



are also plotted in (b). (c) Evolution of distribution function. 
Color scale denotes number of individuals from to 60. 
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FIG. 4: Response for the disturbances added at t*. The 
evolution of distribution function is displayed, (a) Insensitiv- 
ity for small disturbances, (b) Adaptability for large distur- 
bances, (c) Invalid case against large disturbances. 



B. Noiseless case 

We now show the results of numerical simulations. Fol- 
lowing conditions are used in this and the next subsec- 
tions: N 100, a 3.5, p 20.0, tj = 1, ao = 0.1, 
ai = 0.3, a2 = 0.5, 0^3 = 0.2, si = —0.3, S2 = and 
ss = 0.3. Typical snapshots and the evolution of the sys- 
tem are displayed in Fig|2lfor the noiseless case {u = 0). 
As an initial condition, the internal variable Uj{0) of each 
element is randomly distributed and the stock materials 
Wj{0) are all set be zero (FiglH^a)). Each element feels 
the potential U (t) and gathers toward the valley for each 
state. Simultaneously, the elements which staying at the 
j-th valley product the corresponding stock materials Wj 



and change effective depth of the potential. Designed 
ratio (3:5:2) is approximately achieved and the shape of 
the potential becomes almost flat at t = 100 (Fig|3{b)). 
The evolution of distribution function shows the designed 
ratio is almost established in early stage (Fig|3Jc)). 

The response of the system for various type of distur- 
bance is shown in Fig^ First, a few elements in state 3 
are moved to state 2 at as a small disturbance. The 
system is insensitive and reconstruction of distribution 
does not occur (FigEJ^a)). Next, large disturbance is 
given, i.e., all the elements in state 1 are moved to state 
2. In this case, the reconstruction process occurs and the 
ratio between the states are regulated again (Fig^b)). 
Finally, an inadequate example is presented. When all 
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FIG. 5: Typical evolution in noisy case, (a) Noise is added 
from t*. (b) Evolution of randomly selected elements in (a). 



the elements in state 3 are moved to state 2, undesired 
two clusters is realized (Figmc)). This is because the sys- 
tem is trapped in local minimum of the potential. This 
fault can be avoided by the noise term as shown in the 
next subsection. 
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FIG. 6: Recovery from the invalid case shown in Fig^Ic) 
Noise — 4) is added from t*. 




FIG. 7: Ratio design is changed from 3:5:2 to 1:5:4 from t* 
Noise {v — 4) are added from t = 0. 



C. Noisy case 

We now consider noisy case, v ^ {). The noise in- 
troduced here may have physical, chemical or biological 
origin, like the thermal noise of the system, the fluctua- 
tion of distribution of the chemical substances, diversity 
of the age of the individuals . Added noise kicks individ- 
uals away from the potential valley and escape from the 
local minima is expected. Fig |H1 displays a typical evolu- 
tion of the system in noisy case {v = A) . The distribution 
becomes broad but the ratio does not change (FigOa)). 
Fig|5Jb) shows the evolution of randomly selected ele- 
ments as samples of behavior of individuals. During the 
time displayed, most of the elements stay in the same 
state. A few elements, however, change their internal 
state by the noise. The larger v becomes, the more fre- 
quent hopping process between the states occurs. With 
too strong noise, the distribution function collapses into 
single peak. On the other hand, the system hardly es- 
capes from the local minima if the noise is too weak. In 
this sense there is an adequate noise amplitude. Actu- 
ally, noise added from t* set the trapped system free from 
the local minimum as shown in FigEl where the initial 
conditions are set to be trapped case in Fig^c). 



As another example of the adaptability of the system, 
the environment parameter is changed. The shape of the 
potential characterized by (ai,a2,<^3) is changed from 
(0.3,0.5,0.2) to (0.1,0.5,0.4) at T (FiglH). The system 
adapts to the new environment quickly. 

How about the disturbance of changing total number 
N ? We demonstrate an open system including "birth" 
and "death" process as a caricature of "polyethism" ob- 
served in colony of social insect. As a birth process of 
individuals, a new element enters into the system with 
constant internal value u^ew = 0-7 for every five time 
units. As a death process, elements which exceed the 
lower threshold Udeath = —0.31 are removed from the 
system. FiglHl shows the typical evolution of distribution 
function and population of the whole. Both the propor- 
tion between the states and the total number are main- 
tained with fluctuations. 



III. PHEROMONE TRAILING 

So far, various studies have been made on the collec- 
tive behavior of social insect s [l^ll^ll^ . Especially, the 
trail formation of ants (1^ EarEaTMlIl Q have 
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FIG. 8: Example of "open" system with "birth" and "death" 
process. The time series of distribution function (a) and total 
number of elements (b) are displayed. 



called a wide interest as a remarkably synergetic behav- 
ior by presumably not highly intelligent individuals. In 
this section, setting a simple set of rules of which sev- 
eral ingredients are following previous studies 0, , we 
investigate the relation between i)the schedules of feed- 
ing which represent the unsteady natural environment, 
ii)emerging patterns of trails which connect between a 
nest with food resources, and iii)the foraging efficiency 
as a group. Our final purpose is to connect the pattern 
formation of the trail to the group strategy of foraging 
ants under an unsteadily varying environment. 



A. Model 

The model consists of a colony of N movable agents 
(we call them 'ants') which are situated on a 2D hon- 
eycomb lattice with periodic boundary condition. Each 
of these ants is located at one site in the lattice and is 
heading to one of the six nearest sites (Fig. 9) to which, 
the right-hand neighbor or the left-hand, she will move 
according to the dynamical rule described below. The 
task for ants is foraging, that is, finding food and carry- 
ing it back to a nest and the feeding sites are located at 
two corners of a regular triangle with edge-length L and 




possible sites 
in the next step 

site in the 
previous step 



FIG. 9: Each ant is, at each time step, heading one of the 
six nearest sites. This heading direction corresponds to the 
moving direction in the previous time step. The possible 
moving directions in each time step are: the forward (=fac- 
ing) direction and its neighboring (right and left) directions. 
To choose one of them, a stochastic rule with the weight, 
Pa = exp{—A^/T)/Z is adopted where Z is the normaliza- 
tion factor. Here a indicates forward, right, or left direction 
relative to the heading direction, and A"" = p'^—ppiyi^ t) is the 
gradients of pheromone density to the candidate directions 
where (3 is the index to indicate the recruit pheromone (for 
mode-II ants) or the foot pheromone(for mode-Ill ants). Note 
the above weight is not applicable for ants in mode-I which is 
making random walk independent of pheromone field. 



the nest as the third corner (fig. 10). The amount and 
the schedule of feeding are controlled through a pair of 
control parameters. T, M as explained in Fig. 10. 

As the mean for the communication among ants, two 
types of pheromones are secreted and perceived by indi- 
vidual ants according to their temporal modes which con- 
sist of: i)the semi-random walk mode (mode-I), ii)the ex- 
ploring mode (mode-II) and iii)the homing mode (mode- 
III) (Fig. 11). The kind of pheromones considered in 
the present model consists of: i) recruit pheromone and 
ii)foot pheromone. In addition these pheromones evapo- 
rate/decay in a constant rate, also diffuse to the nearest 
sites. 



B. Simulation 



All the simulations are performed on a 150 x 150 hon- 
eycomb lattice with periodic boundary condition. At the 
initial time step, N=500 ants are located at the nest, then 
simultaneously released from the nest. Note the unit time 
in our simulation corresponds to one Monte Cairo step 
within which N randomly chosen ants will sequentially 
take one individual step. Under a proper combination of 
several fixed parameters, e.g., the evaporation/diffusion 
rate of pheromones, and varying a pair of control parame- 
ters {M, T}, the formation process of trail patterns and 
the subsequent efficient foraging are investigated. We 
are particularly interested in the relation between i)the 
geometries of trails and ii)the accompanying foraging ef- 
ficiency. 
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FIG. 10: The method of feeding. Feeding sites are located at 
two corners of an equilateral triangle, the 3rd corner of which 
is the nest site. The amount M of food is supplied, by turn, 
from alternative feeding sites at every feeding interval T. 



(mode-IE) 



(mode-I) 

semi-random walk mode 

secretion of foot pheromone 



exploring mode 

secretion of foot pheromone 
following recruit pheromone 




nest 



(mode-ni) 
homing mode 

secretion of recruit pheromone 
following recruit pheromone 



FIG. 11: Three modes of ants. In mode-I ants make a semi- 
random walk secreting foot pheromone, in mode-II ants ex- 
plore new food following recruit pheromone in which mode 
they secrete foot pheromone, in mode-Ill ant try to come 
back to nest following foot pheromone while secreting re- 
cruit pheromone. The mode switch from mode-I to mode-II 
takes place if an ant reaches a site where the density of re- 
cruit pheromone exceeds a threshold value, whereas a mode- 
Il/mode-III ant turns into mode-III/mode-I on arriving at 
food/nest. 



C. Typical trail shapes 

Once a certain shape of trail is established it is rec- 
ognized as one of following three types i) V-shaped trail, 
ii) Y-shaped trail or iii)/-shaped trail as shown in Fig. 12. 
Notice that these shapes are recognizable and definable 
after averaging the density field of ants over a period 
much longer than the feeding period T. The condition, 
in the M — T space, for the appearance of each of these 
trails are shown in Fig. 13 which indicates that the most 
relevant quantity for the trail geometry seems to be M/T. 
As M/T increases, the shape of trail varies roughly from 
the V-shape to the Y-shape and finally the /-shaped trail 
is obtained. 




• : Food 
X :Nest 

FIG. 12: Typical shapes of trails obtained in the simulation, 
i) V-shaped trail, ii) Y-shaped trail, iii) /-shaped trail. Here 
dark shading in each figure means the density of the ants 
averaged over a time significantly longer than T, the interval 
between successive feeding events. 
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FIG. 13: The relation between emergent trail patterns and 
the combination of feeding schedule parameters, {M,T}. The 
dark shading in each figure means the density of ants averaged 
over a time significantly longer than T and over the ensemble 
of 5 simulations for the V-shaped and the Y-shaped trails. For 
the /-shaped trail the ensemble average is not taken because 
of their intrinsic asymmetry. 



D. Foraging efficiency 

To investigate the relationship between the efficiency 
of foraging and trail geometry, we introduce averaged 
foraging efficiency Eav which is the total amount of food 
carried into the nest per unit time. Hence, the relation 
between M/T and Eav is measured under various feed- 
ing rate M/T (here T is fixed at T = 600). In figure 14, 
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FIG. 14: The relation between the supphed amount of food 
per unit time M/T and the long-term averaged foraging ef- 
ficiency Eav There are five characteristic zones A,A\B,C 
and D. In zone A, the V-shaped trail appears, while in zones 
A' , B the Y-shaped trail appears, and in zones C and D the 
/-shaped trail is seen. Up to the zone A' the perfect foraging 
is attained whereas in zone D the amount of supplied food 
exceeds the carrying capacity of ants thus Eav reaches a sat- 
uration value a. 



a characteristic relation between M/T and Eav is seen 
until Eav reaches an saturation value. In this relation, 
considering the emergent geometry of trails, roughly five 
characteristic zones are recognized as indicated by sym- 
bols A^ A' ^B^C and D. In zone A with the emergence 
of the V-shaped trail, an almost linear ^ — Eav rela- 
tion is seen, and its slope is 1 which corresponds to the 
perfect foraging , namely, ants collect all the supplied 
food at both feeding sites along the V-shaped trail. On 
the other side, in zone D in Fig. 14, ^ — Eav relation 
reaches a plateau, where the averaged foraging efficiency 
Eav amounts to the maximum possible foraging efficiency 
Emax realized by the present number of ants. The last 
situation occurs when the amount of supplied food at 
each feeding site exceeds (or is equal to) the carrying 
capacity of ants, i.e., M/2T > Emax is satisfied where 
the equivalence holds at the left edge of zone D. These 
two zones A and D are characterized by simple forag- 
ing strategies with a straightly extending trail(a pair of 
straight trails) for the /-shaped trail (for the V-shaped 
trail). 

In the case of intermediate amount of food between 
zone A and D, simple foraging tactics by use of straight 
trail(s) fails to give an optimal foraging. For example, 
as M increases over a critical value Mci in zone A' ^ if 
ants stick to the above mentioned simple strategy with 
the V-shaped trail (a pair of straight trails), the sum of 
consuming periods for i) finding new food at alternative 
feeding site and that for ii) carrying away all the food at 
the feeding site, exceeds the feeding interval T. It means 
a part of food remains to be unexhausted even at the end 
of each feeding interval T. However, in our simulation 
at zone A' ^ instead of sticking to such a simple strat- 
egy with the V-shaped trail, ants invent Y-shaped trail 
which, by continuous use of the steadily maintained trunk 
trail near the nest as a common path and leveraging the 



junction of Y as the 'base camp' for the exploration of 
new food, serves the system with a shorter finding time 
of new food at the alternative feeding site. As a result, 
the perfect foraging is kept realized over M = Mci up to 
M = Mc2 in Fig. 14*24]. Certainly, along the Y-shaped 
trail, the distance ants need to walk from the nest to food 
is longer than that along a V-shaped trail. Therefore, 
once a trail is established until all the food at the corre- 
sponding feeding site is exhausted, the temporal foraging 
efficiency along the Y-shaped trail falls below that along 
the V-shaped trail. In this sense, the transition from 
the V-shaped trail to the Y-shaped trail is a trade-off 
strategy within a restricted food environment 24J. This 
perfect foraging of both feeding sites of food is kept main- 
tained before a certain value of M(= Mc2 > Md) over 
which no tactics can realize the perfect foraging. 

In figure 15 we show the theoretically estimated forag- 
ing efficiencies realized by three types of trails as the 
functions of M/T . The figure indicates that the most ef- 
ficient trail pattern varies, as M increases, from V( or Y) 
to Y and finally to /. which estimation, at least qualita- 
tively, corresponds to the simulation outcomes shown in 
Fig. 14. It is remarkable that, even in the imperfect forag- 
ing regime over M = Mc2, through the proper change of 
trail geometry, ants in the simulation keep making the al- 
most optimal foraging close to the theoretically estimated 
one in Fig. 15 under respective feeding conditions. As the 
extreme case in zone-D, the /-shaped trail is built which 
is optimized in the sense that food is sufficient not to run 
out even when ants concentrate on one feeding site and 
in such a case straight trail is most efficient because of 
the minimal path length between the nest and food. Al- 
though the setup of the present simulation is very simple, 
similar arguments to the present are considered to hold in 
more complex and more realistic situations like the trail 
formation under a larger number of feeding sites(Fig.l6), 
in which situations the presumably 'unintelligent' ants 
are guessed to make globally efficient behavior only via 
local information carried by pheromones. 

Still, the basic question; why such an almost optimized 
foraging is realized only with the pheromone trailing tac- 
tics, is left unsolved. However it is strongly expected 
that the further exploration of the pheromone trailing 
will serve us with some intrinsic insights on the reason 
why group behavior of presumably not intellectual ele- 
ments is kept realized in the nature. Also, the applica- 
tion of the pheromone trailing tactics to a wider range 
of systems, e.g., interacting multi-robot system, is one of 
intriguing issue as described in the next section. 



IV. ROBOT EXPERIMENTS 

In previous section, we discuss the mathematical model 
of ant-like system, but we are also trying to develop the 



8 



T=600 



M/T 



FIG. 15: Theoretically estimated relations between M/T and 
averaged foraging efficiency Eav for the V-shaped trail with 
the real line, for the Y-shaped trail with the dotted line, and 
for the /-shaped trail with the broken line, respectively. Ac- 
cording to these relations, the most efficient shape of trail 
varies, as M/T increases, from V(or Y) to Y. Finally /-shape 
is seen most efficient. This way of shape shifting corresponds 
to the emergent trail shapes in our simulation with T=600. 
Note this is the output of a crude theoretical approximation, 
so it contains qualitative deviations compared to those ob- 
tained in the simulation shown in Fig. 14. 
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FIG. 16: Typical trail pattern in the case of four feeding sites 



robots system towards the realization of ant-inspired co- 
operative multi-robot system. 

In engineering field, research in the area of multi-robot 
systems has been very active, and many researchers and 
engineers are currently studying the behavior of robots 
inspired by ant colonies 13j. In most cases, physical me- 
dia such as light, sound, radio wave are used for the com- 
munication between the robots, and some efforts have 
also been made for realizing chemical communication 
robot system 14J. Unfortunately, it is not easy to get 
proper chemical sensors at this stage. Moreover, chemi- 
cal materials, especially gas, are invisible and it is quite 
difficult to observe how they spread and affect robots' 
behaviors. 

In our previous work, we propose "Virtual Dynamic 
Environment for Autonomous Robots (V-DEAR)" for 
real robot experiment 15], in which the pheromones are 
virtually expressed by light information. 



A. Experimental apparatus 

In V-DEAR system, pheromones are replaced with 
computer graphics projected on the floor. Robots decide 
their actions following the color information of the pro- 
jected CG using color sensors. As virtual pheromones are 
represented as CG, we can avoid the problems described 
above. In addition, we can easily control the rate of dif- 
fusion, evaporation, diversity, etc. of the virtual chemical 
materials. 

Fig. [TTT a) shows the schematic of this system. This 
is composed of LC projector to project the CG and the 
CCD camera to trace the position of the robots in the 
field. The robot moving on the field has sensors on the 
top to detect the color and brightness of the field(fig. 
[TTT b)). and determines its behavior. 




FIG. 17: Virtual Dynamic Environment for Autonomous 
Robots (V-DEAR). (a) Schematic of the system, (b) Color sen- 
sors on the top. 

Combining the position information of the robots ac- 
quired from CCD camera and the projected CG by pro- 
jector, we can realize the dynamic interaction between 
the environment and robots. 



B. Task allocation 

Using this system, we made experiments of task alloca- 
tion. In this experiment, each robot can take two states, 
staying and foraging, and decide one of the two states 
based on the total amount of food stock in the nest. The 
stock is consumed constantly depending on the number 
of staying robots in the nest, and it is added when the 
foraging robot returns to the nest. 

Fig. IT^ a) and (b) show the experimental field and 
the schematic, respectively. A half of the field is "the 
nest" and the robots usually move around in this area. 
The other half is "the work space" and the robots mov- 
ing around this area carry back food. V-DEAR system 
monitors the total amount of food stock and expresses 
it by the brightness of the nest. Each robot detects it 
and decides their state autonomously. The robots were 
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adjusted to show nearly 1/3 working ratio in case of six 
robots system (Fig. Efc)). We reduced the number of 
robots by half under the same condition, but you can see 
the rate is kept 1/3 approximately (Fig. [TsT d)). 
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FIG. 18: (a) Experimental field for task allocation, 
(b) Schematic of the experiment. (c)The number of robots 
in each state(six robots). (d)The number of robots in each 
state(three robots). 



C. Pheromone Trailing 

We also made foraging experiments by multi-robot sys- 
tem. Fig. ITUT a) shows the basic behavior of the foraging 
robots. On discovering a food, the robot turns on a LED 
on the top and moves towards the nest. The V-DEAR 
system detects the LED and projects a CG pheromone 
trail during the LED is turned on. When the robot ar- 
rives at the nest, it turns off the LED and changes into 
the searching state. If it finds the pheromone trail, it 
follows the trail. 

Fig. ITUr bVc) show the snapshot of the foraging by 
two robots in case of high evaporation rate and low 
evaporation rate, respectively, and the drawings at the 
right endpoint are the trajectories of the robots. In case 
that the evaporation rate of the pheromone is high, the 
pheromone trail hardly remains. However, the evapora- 
tion rate is low, a stable trail is formed between the food 
point and their nest. 
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FIG. 19: (a)Basic behavior of the foraging robot. On discov- 
ering food, the robot lays chemical trail while returning to 
nest. The robot detecting the pheromone follows the trail. 
(b)In case that the rate of the evaporation is high. (c)In case 
that the rate of the evaporation is low. 



of the colony through the stock materials Wj{t) and de- 
termines its behavior using its internal variables Ui{t) 
and potential U{t). The proportion regulation against 
several types of disturbances was exhibited qualitatively. 
As for general task rule like more complex ones or espe- 
cially asymmetric ones, further improvement seems to be 
required. 

In §3, foraging behavior of ants was simulated, in which 
ants rely only on the local information of two kinds of 
pheromones and move according to a set of primitive 
rules. However, they fulfill a highly efficient foraging 
flexibly adapting to various feeding conditions. Although 
present simulations are restricted to rather simple cases, 
same kind of argument is expected to be available for 
more complex and more realistic situations. 

In §4, small scale experiments of task allocation process 
and trail formation process were examined using multi- 
robot and V-DEAR system which simulates pheromones 
by light information. 

The combination of the above investigations would un- 
cover the fundamental mechanism for the formation, the 
maintenance and the evolution of various types of col- 
lective dynamics exhibited by a wide class of active ele- 
ments. 



V. SUMMARY 

Here in this article, we have investigated two topics 
among intriguing behaviors exhibited by biological mass. 
In §2, we proposed a model for multi(three)-tasks alloca- 
tion in which individuals feel the situation of the whole 
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